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Abstract 

High-mass stars are major players in the chemical and dynamical evolution of galaxies. 
Open clusters and associations represent the natural laboratories to study their evolution. 
In this review, I will present a personal selection of current research topics that highlight the 
use of open clusters to constrain different properties of high-mass stars, such as the possible 
existence of an upper limit for the mass of a star, the evolutionary stage of blue supergiants 
or the characterisation of supernova progenitors. 



1 Introduction 



The study of high-mass stars (also, though not quite correctly, known as massive stars) is 
intimately linked to, and difficult to disentangle from, the study of young open clusters. This 
is not only because most high-mass stars are found within young open clusters and associ- 
ations, but also because clusters are the natural laboratories for investigating the evolution 
of high-mass stars. At the most basic level, determining the fundamental parameters of a 

Since very few 



massive star requires knowledge of its distance (e.g., Herrero et al 



1992) 



high-mass stars have accurate parallaxes (Mai'z Apellaniz et al. , 2008), membership in a clus- 
ter or association is required for calibration. In view of such a close connection, this short 
review can only cover a very small fraction of the many topics subject to active study that 
I could have considered for inclusion. I will concentrate on a few issues of current research 
that highlight the importance of open clusters as laboratories for understanding massive star 
evolution. The selection of these topics is undoubtedly biased by my personal preferences, 
but is representative of the areas currently generating the highest interest among researchers 
in the field. 



2 



Massive stars in clusters 



1.1 Definition of high-mass stars 



There are several possible definitions for high-mass stars, all of them indirect. We can define 
high-mass stars as those initiating C burning non-explosively in their cores. Modern theoreti- 
cal models including overshooting indicate that this will happen for M* ^ 8 Mq (e.g., Eldridge 



Sz Tout , 2004 ) . We can also define high- mass stars as those ending up their lives in supernova 



explosions. These two definitions are almost identical, though the latter implies a slightly 



higher mass (see, e.g., Poelarends et al. , 2008, for a description of the physics involved). The 
limit for a star to produce a supernova explosion has recently been set from observations of 

and see also Sect. [s]). A third 



supernova progenitors at ^ S.S^J gM© (Smartt et al. 



2009 



possible definition of high-mass stars is those stars with self-initiating radiation-driven winds 
(Kudritzki &: Puis, 2000). Radiative winds become detectable for main-sequence stars with 



spectral type earlier than B2V, which again roughly corresponds to ^ 8Mq. 

Observationally, high-mass stars comprise the OB stars (approximately, 02~B2 V, 02- 
B5 III and 02-B9 1; |Reed|2009t [Walborn fc Fitzpatrick|1990l ) and some later type supergiants 
(the most luminous supergiants of spectral types A, F, G, and K, and the M-type supergiants). 



For an observational review, see Massey (2003). 



1.2 Formation of high-mass stars 

The formation of high-mass stars is a major research topic in modern astrophysics, and will 



not be discussed here. See Zinnecker & Yorke (2007) for a recent review. High-mass stars 



may have an effect on the formation of low-mass stars, and thus on the Initial Mass Function 
(IMF). In addition, high-mass stars may play a role in triggering further star formation. This 



issue has been controversial since the emergence of the classical theory ( Elmegreen Sz Lada 



1977), because a causal relationship is difficult to assess, even if indications of sequential 



formation are widespread and strong indirect evidence for triggering has been found (e.g.. 



Walborn, 2002 Zavagno et al. , 2005). I will not discuss star formation here, but simply note 



that sequential star formation within a cluster or association may give rise to populations 
with different ages, which might be difficult to disentangle, meaning that accurate ageing is 
not always possible for young open clusters. 



2 The most massive stars 



The possible existence of an upper limit to the IMF has been a hotly debated issue in recent 
times. Several stars believed to be extremely luminous, and thus extremely massive, have 
been resolved into two or more components (e.g., Mafz Apellaniz et al. , 2007). Based on an 



analysis of the stellar population in the compact young Arches open cluster, Figer (2005) 
concluded that there was clear evidence for an upper limit to stellar masses around 150 M©. 
This result is dependent on several factors, such as the extinction law to the Galactic Centre, 



but has also been reproduced for the Large Magellanic Cloud (Koen, 2006) 



Recent modelling advances have revealed that hydrogen-rich Wolf-Rayet (WR) stars of 
the nitrogen sequence (WN) are really main-sequence objects with heavy mass loss. These 



are candidates to be the most massive stars in their hydrogen core burning phase (Langer et 
al. , 1994 de Koter et al. , 1997 Crowther Sz Dessart 1998). The fact that they are commonly 
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found heavily packed at the core of the densest young open clusters supports this view. 
Searches for very massive stars have thus concentrated on possible binaries amongst this 
population, as close binaries may provide accurate masses via dynamical studies. 

One of the best examples is found in the obscured starburst cluster Westerlund 2. 
Located in the Sagittarius Arm (likely at ~ 8 kpc, though distances as short as 3 kpc have 
been given), this cluster contains at least 14 early O or WR stars (Rauw et al. , 2007). Its 
total mass is also open to discussion, but may approach Md = 10^ Mq ( Ascenso et al. , 2007), 



a result supported by the recent discovery of two very massive stars likely ejected from the 
cluster core (Roman-Lopes et al. , 2011 ). One of the most luminous members of Westerlund 2, 
the WN6h star WR20a, was identified as a massive binary by Rauw et al. (2004) and as an 
eclipsing binary by Bonanos et al. (2004). Combination of the radial velocity curve and the 
lightcurve provides an accurate solution with Porb = 3.69 it 0.01 d. The components have 
masses Mi = 83 it 5 Mq and M2 = 82 it 5 Mq . They have the same spectral type and the 
same mass, which is the highest stellar mass measured with high accuracy. 

Even more massive than Westerlund 2, the cluster at the core of the giant Hll region 
NGC 3603, also located in the Sagittarius arm (at ~ 7 kpc), is extraordinarily compact. 
It contains ~ 35 early O or WR stars (Moffat et al. , 1994) and likely has a mass Aid = 



1.5 X 10^ Mq (Rochau et al. , 2010, and references therein). Its central concentration contains 
3 WN6ha stars, the brightest of which is cin eclipsing binsiry with. .Forb — 

3.7724 d. The 

orbit has been solved using K-hand spectroscopy (Schnurr et al. , 2008). The radial velocity 
curve of the secondary has large uncertainties and therefore the masses are not very tightly 
constrained. The solution implies Mi = 116 it 31 Mq and M2 = 89 ± 16 Mq. Component 1 
is thus the most massive star with a dynamically determined mass (Schnurr et al. , 2008). 

Outside the Milky Way, the 30 Dor complex in the Large Magellanic Cloud is likely 
the most massive starburst in the Local Group. Its nuclear cluster, Radcliffe 136, is 2.7 pc 
across and contains hundreds of OB stars. Its mass has been estimated to be, at least, 
10^ M0 (Andersen et al. , 2009). Using updated stellar models, Crowther et 



M. 



1 X 



al. (2010) found that the dynamical masses of the WN6ha stars in NGC 3603 could be well 



reproduced by evolutionary tracks. If these models are then extrapolated to the even brighter 
WN5h stars at the core of 30 Dor, they indicate enormous masses. Three of the WN5h stars 



are fitted with models implying current masses > 150 Mq, with Star al having 265+35 
(Crowther et al. , 2010). This would be the most massive star known, well above the proposed 
upper mass limit. Unfortunately, none of the stars at the core of 30 Dor appears to be a 
spectroscopic binary to confirm this determination with a dynamical measurement. 



3 Stellar evolution 

High-mass stars play a decisive role in driving the chemical evolution of galaxies. Heavy mass 
loss through all their life stages and in the final supernova explosion provides an important 
fraction of the heavy elements in the interstellar medium. In order to understand the chemical 
enrichment of the medium, we need to constrain how and when mass is lost by high-mass 
stars, and this means being able to map observed phases on to theoretical evolutionary tracks. 

Open clusters and associations allow us to explore the evolutionary context of massive 



stars in different evolutionary stages (e.g., Walborn, 2010). This process starts with the 
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youngest stellar systems. For example, the Carina nebula provides many of the MK standards 
for the earliest O subtypes (Walborn et al. , 2002 1 . The fact that 03-4 stars are still on the 
main sequence indicates that the complex is very young, while this very youth constrains its 
only evolved star, r] Car, to be extremely massive. 

At a slightly older age, Cyg 0B2 represents an important laboratory for the study of 



0-type stars. Considered a very promising target since its discovery (Johnson Sz Morgan 



1954), it could only be studied in detail after the advent of CCD detectors because of the 



high obscuration, in spite of its relative low distance (see Fig.[T]). [Massey &: Thompson ( 1994) 
identified ~ 60 stars more massive than 15Mq, finding that Cyg 0B2 is very compact for an 
OB association and seems to occupy a more or less unique position somewhat intermediate 
between an open cluster and a normal OB association (cf. Knodlseder 2000). The central 



region contains two cluster- like stellar concentrations (Bica et al. , 2003). They form an 
elongated figure ~ 4' x 10' (~ 2 x 5 pc at the distance of Cyg 0B2) and contain 18 0-type 
stars (assuming that all systems are single; given the observed multiplicity - see below -, the 
actual number is likely to be a factor of two higher). 

After analysis of the 2MASS data in this region suggested a much richer population of 



early type stars (Knodlseder 2000), several recent studies (e.g., Comeron et al. , 2002 Hanson 



2003) have shown that the association contains > 50 (and likely ~ 70-80) 0-type stars (e.g.. 



Kiminki et al. , 2007 Negueruela et al. , 2008 ). The main sequence is very well defined, showing 



a clear turn-off at 06 V. A number of giants and supergiants follow the 2.5 Myr isochrone 
from this turnoff. However, there is also a population of 03-05 supergiants, which seems 
incompatible with this age (see Negueruela et al. 2008, and Fig[T]). With the data available, 
it is impossible to rule out the possibility that these stars represent a younger population. 
However, the similarity of the HR diagram for Cyg 0B2 to theoretical plots in which very 



fast rotation generates blue stragglers (de Mink et al. 2009) is striking 



The presence of blue supergiants with spectral types too early for the main-sequence 
turnoff is very frequent in young open clusters, but these objects have not generally been 
classified as blue stragglers, even if they are located to the left of the cluster isochrone, the 
classical definition of a blue straggler. A paradigmatic case is h Persei (NGC 869), one of 
the most well-studied clusters in the Northern sky, and a rather massive cluster - recent 
mass estimates indicate ~ 4 500Mq in the core region (Currie et al. , 2010). Most recent 
works seem to converge on an age ~ 14 Myr, in good agreement with an apparent main- 
sequence turn-off at BIV. The luminous supergiants in the cluster, however, have spectral 
types too early for this age (Marco Sz Bernabeu 2001), at B2Ia (HD 14143) and B3Ia 



(HD 14134), and evolutionary masses approaching 40 Mq (McErlean et al. , 1999 



et al.[ 2006a). A similar situation is found in many other clusters (Marco et al. 



Crowther 



2007). In 



most cases, there are no observational reasons to suggest that these supergiants are members 
of a younger (second-generation) population. This clearly shows that we are still very far 
from understanding exactly which evolutionary phase blue supergiants represent. 

Because of the huge size of its stellar population, the young cluster Westerlund 1 (Wd 1) 
is currently our best laboratory for studying the evolution of high-mass stars. Located at 

2011) and may approach 10^ Mc 



5 kpc, its mass exceeds 5 x 10 Mq (Gennaro et al. 



0- 

With more than twenty WR stars ( Crowther et al. , 2006b ) and several dozen OB supergiants 
(Negueruela et al. , 2010), Wd 1 provides stringent tests on current theoretical models. The 
population observed is consistent with a single burst of star formation. There may be some 
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Figure 1: Semi-empirical HR diagram for Cyg 0B2, after Negueruela et al. (2008). Note 
that, in the original paper, the {V — K)q term was added twice, because of a mistake with the 
spreadsheet. When this is taken into account, the distance modulus DM = 10.8 proposed 
by Hanson (2003) is no longer favoured. A better fit is now obtained with the classical value 
DM = 11.3, as displayed. Continuous lines are non-rotating isochrones for logt = 6.2, 6.3 
and 6.4 from Schaller et al. (1992). The dashed line is the logt = 6.4 isochrone in the high- 
rotation models (Meynet &: Maeder, 2003). The main-sequence turn-off at 06 V seems to be 
followed by a number of evolved stars tracing the logt = 6.4 isochrone. A number of 03-5 
supergiants, though, fall well to the left of the isochrone. 
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mild blue stragglers, but the blue supergiants in the cluster (including some very luminous 
late-B hypergiants with My ~ —9) fit rather well the theoretical isochrones for ages between 
5 and 6 Myr ( Negueruela et aTj 2010). Moreover, the dynamical determination of the mass 
for the two components of the eclipsing binary Wdl-W13 shows good agreement with the 
expected theoretical masses, with the Bl I component having a mass Ad^, = 35 it 5 Mq (Ritchie 



et al. , 2010). This good agreement is comforting in view of the difficulties in understanding 



blue supergiants in other clusters and provides strong support for the basic assumptions 
underlying current evolutionary tracks (e.g., Meynet & Maeder, 2000). 



4 Binarity 



The large number of high-mass stars present in Wd 1 offers an ideal opportunity to study 
their binary fraction. The population of WR stars offers strong indirect evidence for a very 
high binary fraction ( Crowther et al. , 2006b ) . Targeted observations of the blue supergiant 



population revealed that at least 40% of the stars observed are binaries (Ritchie et al. , 2009). 



Observations of a larger sample, including stars closer to the main sequence, are in progress. 
Searches for binaries in Cyg 0B2 are, to date, more complete. A radial velocity survey 



carried out over several years (Kiminki et al. , 2007 Kobulnicky &; Kiminki, 2011) should 



be able to reveal, except for very unfavourable inclinations, all massive companions and a 
substantial fraction of low-mass companions. The survey has so far detected 20 spectroscopic 
binaries, while 20 other objects show radial velocity variations, and might be binaries. Com- 
parison to theoretical expectations suggests a very high binary fraction, perhaps approaching 



unity. The masses of the companions seem to follow a flat distribution ( Kiminki et al. 2009 



Kobulnicky & Kiminki, 2011). 



The binary fraction and companion mass distribution do not seem to be universal, 
though. Variability in both observables is high among the few clusters with dedicated studies. 
In NGC 6231, the complete sample of 0-type stars (16 objects) presents /bin > 0.63, with 
companions having consistently M^.. > 5 Mq (Sana et al. , 2008). In contrast, the smaller 



cluster NGC 2244, with 6 0-type stars, has a much smaller binary fraction, /bin > 0.17, 
while in IC 1805 (10 0-type stars) /bin > 0.20, with the possibility of massive companions 



almost ruled out (Mahy et al., 2009 de Becker et al. 2006). 



Meanwhile, spectroscopic monitoring of a large sample of O and WN stars has revealed a 



very high fraction of spectroscopic binaries (Barba et al. , 2010). Interferometric observations. 



which sample a different range of orbital sizes, have also demonstrated a very high degree 
of multiplicity amongst high-mass stars. 0-type stars in clusters display a much higher 



binary fraction than field O stars (Mason et al. , 2009). This can be naturally explained if 



many (most) field stars have been ejected from open clusters (in accordance with theoretical 
predictions presented in several contributions to these proceedings). 



5 Cool massive stars 

After the blue supergiant phase, the evolutionary tracks of massive stars become more un- 
certain. An important fraction of high-mass stars (perhaps all with ^ 25 Mq) become 
WR stars. Before reaching this hydrogen-depleted phase, the stars must lose their envelopes. 
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but the evolutionary phase when this mass loss takes place is unclear. Wd 1 contains a high 
number of yellow hypergiants (with spectral types A and F), which seem to be subject to 
heavy pulsations (Clark et al. , 2010). Current evolutionary tracks would predict these objects 



to become Luminous Blue Variables (LBV) and then reach the WR stage. Nevertheless, the 
cluster also contains a number of red supergiants, which are not predicted to exist at such 



high luminosities (Clark et al. 2010). 



Stars less massive than ~ 25 Mq are expected to end their lives as red supergiants 
(RSGs). Modelling of RSGs is very complex due to a number of factors, among which we 
can cite their huge size, poorly determined molecular opacities and the expected heavy mass 
loss. Because of these difficulties, their fundamental properties are not very well known, 
though important advances have been obtained with new model fits to a sample of RSGs in 
open clusters and associations (Levesque et al. , 2005). Spectroscopic monitoring of RSGs in 



clusters has allowed the detection of irregular radial velocity variations, which may be the 



signature of pulsation, the most likely source of mass loss (Mermilliod et al. , 2008). 



In recent years, RSGs have shown their usefulness as signposts of massive clusters. As 
they are very bright in the infrared, they can be seen through heavy absorption. Several 



clusters rich in RSGs have been discovered towards the inner Galaxy (e.g., Figer et al. , 2006 



Clark et al. 2009a). Determination of their chemical composition has provided extremely 
valuable clues to the Galactic chemical distribution (Davies et al. , 2009). These studies are 



allowing the discovery of an increasing number of massive young open clusters in the Galaxy 
(see Portegies Zwart et al. , 2010, for a review). Unfortunately, the actual masses of these 
obscured clusters are not known, as we can only see the RSGs and not the associated main 
sequence (see Fig. 2|) . The most extreme case is the open cluster Stephenson 2, which contains 
at least 26 RSGs (Davies et al. , 2007), and might have a mass approaching 10^ Mq, if current 
predictions for the duration of the RSG phase are correct. Unfortunately, this is not certain, 
as the case of the open cluster NGC 7419 (Marco & Negueruela, 2011, and Fig. [2]) highlights. 

One important test for stellar evolution models is the ratio of blue to red supergiants 
in a population. This ratio is very sensitive to input physics, such as treatment of mass 
loss, convection and mixing processes. Unfortunately, in most Galactic clusters, the number 
of supergiants is so small that the ratio is not statistically significant. Studies so far have 
considered average values over a number of open clusters, finding an apparent dependence 
on Galactocentric radius. When observations of the Magellanic Clouds are also included, 
there is strong evidence for an increase in the blue to red supergiant ratio with increasing 
metallicity (Eggenberger et al. , 2002). This is exactly the opposite behaviour to what current 



evolutionary tracks would predict (Meynet et al. 2011) 



However, simply counting all red and blue supergiants may provide a distorted picture, 
as not all "supergiants" sample the same population. As an example, the open cluster 
NGC 6649, located at ~ 2 kpc in the Sagittarius arm, has an age ~ 50 Myr ( Turner| 1981) 
and contains at least three low-luminosity RSGs (with spectral types around Kl lb) and 
one supergiant Cepheid variable. According to modern evolutionary tracks (e.g., Marigo et 
al. , 2008), the progenitors of these supergiants had initial masses ~ 8Mq. Stars of such 



low masses should never appear as blue supergiants (according to the same tracks), and so 
including them in counts to determine the red to blue ratio may be misleading. 

Nevertheless, these low-mass RSGs are very interesting objects in themselves, as they 
must be the progenitors of most Type II supernovae. This is a simple consequence of the shape 
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Figure 2: Left: 2MASS image of the open cluster NGC 7419. This obscured cluster in 
the Perseus Arm contains 5 red supergiants for no blue supergiants and presents the highest 
fraction of Be stars among Galactic clusters (Marco &; Negueruela 2011). Current theoretical 
models predict that this cluster should have a mass approaching 10^ Mq to contain 5 RSGs 



( [Clark et al.| [2009bj), but its mass is unhkely to exceed 3 x lO'^ Mq. Right: 2MASS image 
of RSGCl, a cluster of red supergiants at d ~ 6 kpc, in the Inner Galaxy (Figer et al. 
20061). This cluster has an age ~ 10 Myr and may have > 3 x 10"^ Mq (iDavies et al.l |2008|). 



Comparison of this image to the one in the left panel highlights the brightness of RSGs in 
the near infrared and their role as signposts of recent star formation. The population of blue 
stars easily seen in NGC 7419 fades out in RSGCl due to the higher distance, much higher 
extinction, and confusion. 
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of the IMF, as there must be many more stars with M^, ~ 8Mq than with ~ 20 Mq. As 
Type II supernovae are very important contributors to the dynamical and chemical evolution 
of galaxies, determining the lower mass for a supernova explosion becomes a fundamental 
issue to understand the history of the Universe. In recent years, an important observational 
effort has been dedicated to identifying supernova progenitors and deriving their properties 
(see Smartt, 2009). In some cases, it has been possible to identify the actual star in pre- 
explosion images. Sometimes, the supernova seems to have taken place within an open 
cluster or association and properties can be inferred from the parent population. 

Based on these o bservations, the min imum mass for an exploding star has been esti- 
mated at S.S^i gMo (Smartt et al., 2009). Efforts are also being dedicated to determining 
whether there is a direct relation between the initial mass of a star and the class of super- 
nova explosion it will undergo. Of course, the final fate of a massive star depends on several 
factors, such as mass loss rates during the different evolutionary phases, initial rotational 



velocity and its evolution, and, perhaps most decisively, binary interaction (e.g., Meynet et 



al. , 2011). In spite of this, attempts have been made to the develop "typical" scenarios (e.g.. 



Smith et al., 2011) 



Another interesting fact is the apparent correlation between the fraction of Be stars 
among B-type stars and the number of RSGs. Both seem to increase in a similar way with 
decreasing metallicity (Meynet et al. , 2007), suggesting a possible evolutionary connection 
with fast rotation. In this sense, the open cluster NGC 7419 (Fig. [2]), which has the highest 
fraction of Be stars amongst Milky Way clusters, presents a blue to red supergiant ratio of 0/5 
(Marco &: Negueruela, 2011), atypical for Milky Way clusters. In contrast, NGC 663, which 
also has a very high Be-star fraction and is located at approximately the same Galactocentric 
distance, presents a ratio 5/0 for the core region (6/2 when the halo is also considered), 
showing that correlations are not always direct or simply that statistical fluctuations may 
dominate the observables even in moderately massive clusters. 



6 Outlook 

In the near future, Gaia is bound to play a decisive role in furthering our knowledge of high- 
mass stars. Accurate parallaxes will result in much better distances to open clusters, implying 
improved luminosity determinations. In addition, confirmation of membership for individual 
peculiar objects will give them an evolutionary context, providing strong constraints for 
evolutionary tracks. Gaia may also contribute strongly to the study of binarity amongst 
high-mass stars, as it will obtain accurate lightcurves for huge samples and may detect orbital 
motions in wide binary systems. 

Of course, Gaia cannot solve every problem. We will need other data and also new 
theoretical developments. High-resolution spectroscopy coupled with accurate stellar atmo- 
sphere models will be needed to determine stellar parameters beyond luminosity (see, e.g., 
Simon-Diaz et al., in these proceedings). Improved evolutionary tracks will result from deeper 
interaction between theory and observations. Radial velocity curve solutions for spectroscopic 
binaries will sample the intermediate range of binary separations. 

But, above all, we have to be conscious of one main limitation: Gaia will not be able 
to sample the inner Galaxy, where high extinction and crowding will impede the acquisi- 
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tion of high-quality data. Almost all the massive young clusters known will be beyond its 
reach. Fortunately, a full complement of new instruments will help us to study star formation 
in the inner Galaxy. Among them, the first generation of infrared spectrographs with high 
multiplexing capabilities, such as KMOS at the VLT or MIRADAS at GTC, will play a funda- 
mental role in studying obscured clusters and massive star-forming regions. Meanwhile, new 
space telescopes will provide the tools to resolve clusters in the galaxies of the Local Group. 
JWST may be an ideal tool for detecting RSGs in distant galaxies (and hence identifying 
supernova progenitors), while WSO/UV will be used to study low-extinction populations of 
massive hot stars. In all, with this formidable set of new missions and instruments, the next 
few years are likely to see important developments in this field. 
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